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Abstract: With reference to a particle tunneling through two successive 
barriers, it seems to have been generally accepted that the tunneling time 
does not depend on the separation distance between the barriers. This phe- 
nomenon has been called the generalized Hartman effect. In this letter, we 
point out a lack of mathematical rigour in the reasoning by which this ef- 
fect was deduced about ten years ago. A mathematically rigorous treatment 
shows us that the tunneling time does indeed depend on the length of the 
free space between the barriers. 
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Consider the integral 

1 = j G{E)e''^^UE, (1) 

where G{E) and 0{E) are real valued functions of a real variable E, and G{E) 
is assumed to be a sharply peaked function such as a Gaussian centered at a 
mean value. Denote the deviation of G{E) about its mean value Ei as 6E. 
Let 

9{E) = e{Ei) + e'{Ei){E - El) + 0[{E - Eif] (2) 

be a Taylor expansion of 9{E). In the case that we can neglect the third 
term of the right side above for some reason, the integral / can be rewritten 
as 

/ = e^^(^^) j G(E)e*^'(^^)(^-^^)rfE. (3) 

If 

e'{Ei)5E > 1, (4) 



^ Corresponding author. E-mail address: shuichi@edu.u-ryukyu.ac.jp 
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Figure 1: G{E) given by Eq. = 5) 



then the term e*^'(^i)(^ -^i) rapidly oscillates as E varies around Ei, and / 
approximates to zero. In contrast, if 6'{Ei)6E ^ 0, especially if 

e'{E,) = 0, (5) 

then the oscillation is suppressed and / results in an appreciable value. 

This is the essence of the stationary phase method (SPM). Here, we have 
to emphasize that the condition indicated in italics for the function G{E) is 
very important in the reasoning. For instance, consider the function 

G{E) = e-^' cos(rE) (6) 

where r is a real parameter. Fig. [1] shows the graph of G{E) for r = 5. In 
that G{E) firstly takes both positive and negative values and secondly has 
multiple extrema, this function does not satisfy the condition indicated in 
italics. 

Using G(E) given by Eq. (E]), consider the integral 

I{t) = f G{E)e''^dE (7) 



where t is a real parameter. If we apply the SPM to this integral, since 

^tE = t, (8) 

we are led to the conclusion that the value of I(t) is appreciable only for 
t = and is negligible for sufficiently large t. 
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Figure 2: I{t) given by Eq. (r = 5) 

However, since 

G{E) = e-^' (e'"^ + e"'"^) /2, (9) 
I(t) can be expressed as 

/(t) = ^e-^t+r)y4 ^ g-(t-.)V4) . (10) 

The t for which the value of I{t) is appreciable is therefore not t = but 
t = ±r. The graph of I{t) as a function of t is shown in Fig. [2] for r = 5. 

Thus, when we use the SPM for the integral ([T]), we must not forget to 
check whether G{E) is a sharply peaked function such as a Gaussian. 

Now, consider a particle tunneling through two equal successive rectan- 
gular barriers with height Vq and width a: 

V{x) = Vo [e{x)e{a -x) + e{x - L)e{L + a-x)] (11) 

where < a < L, and L — a is the distance between the barriers. About 
ten years ago, Olkhovsky, Recami and Salesi [1] claimed that the tunneling 
time does not depend on the distance between the two barrier^. However, 
we conclude that: 

1. Although the SPM is used in the essential part of their reasoning in 
[TJAits applicability is questionable. The doubt arises in that their use 
of it is essentially the same as the incorrect use exemplified above. 

^ For an arbitrary number of barriers, see [2]. 
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2. When we correctly apply the SPM to the wave packet which describes 
the particle, the tunneling time does in fact depend on the distance 
between the two barriers. 



3. Contrary to our predictions, Longhi, Laporta, Belmonte and Recami 
[3] claim to have confirmed with experimental evidence the theoretical 
predictions obtained in Ref. [T]. If this independence of the barrier 
separation is indeed confirmed, we are faced with a serious contradiction 
between our theory and physical reality. 

First, consider the one-dimensional Schrodinger equation 

1 (P 



2m dx^ 



+ {E-V{x)) 



i^E{x) = 0. (12) 



If we consider the plane wave e*^^' with k = \/2mE as the incident wave to 
the left of the potential, then we have 

^e{x) = T(E)e^^^ {x> L + a) (13) 

as the transmitted wave to the right of the potential, where the coefficient 
T{E) has the form (see p. 52 of Ref. |1]) 

„— ifc(L+a) 

(14) 

(cosh^a + iA_ sinhxa)^e~*'^(^~'^) + (A_|_ sinhxa)^e*'^(^~'^) 

with 

X = ^2m{Vo - E) 

and 

Integrating over a band of stationary states with different energies, we 
construct a spatially localized incident wave packet 

g{E)e'^''-'^'dE, (15) 

for which the transmitted wave has the form 

j g{E)T{E)e'^''-''^UE. (16) 



In the following, g{E) is assumed to be a Gaussian distribution given by 



9iE) 



~{E~Eo)^/25^ 



;i7) 



where Eq is within the interval (0, Vq) and 5 is rather small, such that the 
majority of the energy distribution is contained in (0, Vq). 

The authors of [1] restrict themselves to the case where a is large enough 
(and X iiot too small) that one can assume xa — )■ oo, and they obtain 



where 



A 



[ik — x) 
2kx 



(18) 



(19) 



2kx cos k{L — a) + (x^ — k"^) sin k{L — a) 

After confirming only that A is real, they immediately conclude that we can 
derive that the tunneling time 



M:un 



daig T{E)e 



ik{L+a) 



BE 



d_ 
dE 



arctan 



2kx 



_d_ 

dE 
2m 
kx 



arg 



—Aikx 
{ik - xf 



(20) 



while depending on the energy of the particle, does not depend on L + a 
(being actually independent of both a and L). (See p. 882 of Ref. [Ij.) 

If the term tunneling time means the time when the peak of the trans- 
mitted wave emerges at the point x = L + a, then their conclusion seems to 
be in doubt. The reason is as follows: In the limit XQ- — C)0, the transmitted 
wave given by Eq. flTBl) has the form 



g{E)T{E)e'^''-'^'dE 



g{E)e-^''''A- 
-kx 



-4ikx 



{ik - xy 

[ (^(E)e-2^'^A-i^e^{^(^-^~")+^-^*>dE, 
J k -\- y 



^i{k{x-L-a)-Et} ^j^ 

(21) 



where 



-4ikx 



{ik - x)^ k^ + X^ 



9 = arctan ■ 



2kx 



(22) 
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Figure 3: g{E)e-'^^''A^^ given by Eq. ([23 



// the function 



g{E)e-^'"'A- 



Akx 



(23) 



is of Gaussian type, then using the SPM we can conclude that the position 
of the peak of the transmitted wave is given by 



dk 



de 



— ix - L-a) + — - t = 0. 



dE 



dE 



This imphes that the peak emerges at x = L + a at time 

_ de 

^~dE' 



(24) 



(25) 



that is to say, the equation (1201) is reasonable. 

But the function (!23l) is far from Gaussian: Fig. [3] shows the graph of 
(i23|) with the width a = 5/y2rn% and Eq = 5 = Vq/IO. The scale on 
the horizontal axis indicates the ratio E/Vq. The left-hand figure is for the 
case L — a = 8/ A/2ml^ and the right-hand figure for L — a = 9/^/2mVo. In 
that it has both positive and negative values and exhibits multiple extrema, 
the function fl23|) is similar to the function (|6]). Therefore, as exemplified at 
the beginning of this paper, there is a possibility that r/^n given by Eq. fl20l) 
does not represent the true tunneling tim^ . 

One further question arises in the reasoning given in [1]: Even if the 
function ( l23l) is like a Gaussian centered on a certain point Ei, the tunneling 



This kind of way of application of the SPM is seen frequently in literatures on this 



area. 
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Figure 4: dO/dE = 2m/ kx 

time is not just the term dO/dE of Eq. fl25|) . dO/dE is a function of E, and 
the tunnehng time t^^^ is obtained only after substituting Ei for its variable 
E. Hence, even if the form of the function dO/dE is independent of both a 
and L as shown by Eq. fl2U]) . the tunneling time may depend on either a or 
L in the case that Ei depends on a or L. 

Since the function (l23l) is far from Gaussian, we cannot define E'l as its 
mean energy. But, interpreting things broadly, we might try to define Ei as 
the point at which the function tends to infinity. The value of A is infinite if 
and only if its denominator 

2kx cos k[L — a) + (x^ — /c^) sin k{L — a) 

becomes null. Those null points vary as L — a changes and the variation is 
rather sensitive to the change of L — a as shown in Fig. [3l and moreover 
those points occasionallv come near the endpoints of the interval (0, Vq). As 
we can see from Fig. the gradient of the function dO/dE is very steep in 
the vicinity of those endpoints. Hence, there is a rather high possibility that 
the value of 

^(E,) (26) 

depends on L — a. 

Next, we show what we can obtain from a cautious application of the 
SPM to the transmitted wave (IT^ : We do not use the approximation (IT5]1 . 
but use the original form fll4p . Introducing r and by 

cosh xa + i A_ sinh xa = re"**^, 
^ The scales on the horizontal axis indicate the ratio i?/Vo. 
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1 + sinh xa, = — arctan(A_ tanh^a), (27) 
we have 

„—ik{L+a) 

(cosh xa + 2A_ sinh xa)^e~**^(^~") + (A+ sinh xa)^e*'^(^^") 

(A+ sinh xa)^ 



^—ik{L+a) 



^ ^-ik(L+a) ^ 1 / A+Sinhxa \ ^i{(2„+l)fc{L-a)+2{n+l)<;i} 

The transmitted wave f|T6|) therefore can be written as 

2n 



1)"! ^^^^ I^ A+sinhxa 

^ gj{fc(x-L-a)+(2ri+l)fc(L-a)+2(n+l)</.-£;t} ^29) 



n=0 



that is to say, it is an infinite sum of the wave packets parameterized by 



n = 0, 1, 2, ■ ■ ■ each of which has 



g{E) ( A+sinhxa~ 



2n 



(n = 0,l,2,---) (30) 



as its amphtude. 

Fig. |5] shows the graphs of fl30l) for the cases n = (upper), 10 (middle) 
and 20 (lower), where we set Eq = Vo/2,6 = Vo/10,a = 5/^/2rnVo, and the 
scales on the horizontal axis represent the ratio -E/VoB 

As shown in this figure, as far as a does not exceed some critical valu^, 
each of the amplitudes (130|) is a sharply peaked function centered at a certain 

^ Denote the transmission and the reflection coefficients for one potential barrier of 
width a as To and Rq, respectively. Then we have [4] 

To — (coshxa + iA_ sinh^o)^^, ^o = — jA+ sinhxa(coshxa + iA_ sinh^a)"'^ 

and 

m=9iE)\To\^\Ro\^''. 



^ The transmission for a very thick barrier is over the top of the barrier for the most 
part, essentially not tunneling at all. This fact is stated as one of the conclusions of 
Hartman [5] , and the critical value of a depends on Eq and S determining the distribution 
g{E) defined by Eq. ([TT]) . We should pay attention to this fact when we treat the limit 
a oo. (See [B].) 
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Figure 5: Graphs of Eq. fl30|l 

point El slightly higher than Eq. Therefore, the SPM can be applied to each 
term of the sum (12^ . and the position of the peak corresponding to n is 
given by 

^{Ei){x -L-a) + {2n + l){L- a)^{E,) + 2{n + 1)^{E,) -t = 0. 
Hence, the n-th peak of the transmitted wave emerges at a; = L + a at 



t={2n + 1)(L - a)^m + 2{n + 



(31) 



Since L is not included in (130 p . Ei does not depend on L, and dk/dE{Ei) 
and d(f)/dE{Ei) are hence independent of L. Thus the time t given by Eq. 
fl3T|) is in proportion to L — a for each a and n. 

Define v and r by 



1 dk m 
V ~ dE^ ~ ^f2^i 



(32) 



where (see p. 9 of Ref. ^] 



dd) ma / . r, , , 

'k_ ^ x\ tanh xa 
,X kj xa 



k A / 



1 



cosh xa 



(33) 
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and 



dd) , ^ . 2m 
dE^ k{E,)x{Ei) 



(34) 



The time r is the tunneling time of a particle tunnehng through one barrier, 
and it is independent of the barrier width under the condition, as shown by 
Hartman [5pl . that the barrier is neither very thin nor very thick. 

Eq. ( 13T]) means that the peak of the n-th component of the transmitted 
wave ( 129|) emerges at x = L + a at 



t = {2n + l)'- 



+ 2(n + l)r. 



These times do thus depend on the distance L — a. 
If we consider the behavior of the wave 



(35) 



J g{E)'il^E{x)e-''^'dE 



(36) 



in the regions x < and a < x < L, we can gain a better understanding 
of why the transmitted wave has muhiple peaks: We define the coefficients 
R{E),a{E) and P{E) by 



i>E{x) 

Simple calculations show that 
RiE) 



(37) 



-zA+ sinh xa 



(cosh xa + iA„ sinh xa) e ^^^^ + sinh xoe*^(-^ 
X (cosh xa — ^ A_ sinh xa) e**"^^""^ + (cosh xo- + ^A_ sinh xa) e"*''^^""^ 

Using r, (j) introduced by Eq. flTTj) . we have 

A+sinh^Xa 2i 



/^(i?) = — zA+sinhxo. 



(re-^'^)2e-''=(^-'*) 1 + 



ik(L—a) 



(re-*<^)^ 



X 



^^ik{L-a)+iij> _|_ ^^-ik{L-a)-i(j> 



00 

n=0 



A+ sinh xo 



2?i+l 



^ In his terminology r = (5t3 
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X 



J{2{n+l)k{L-a)+2{n+l)(f)} , i{2nk(L-a)+2n<j)} 



.A+sinhxa 

-t e ^ 

r 

oo 



2n-l 



n=l 



The decomposition fl38l) implies that the reflected part of the wave (|36|) 
is an infinite sum 

j^l'^(^_lY I f A+Smh ^i{-kx+2nk{L^a)+C2n+l)(f>-Et}^j^ 

where the amplitudes 

^(^) and j (n = l,2,---) (39) 

are all Gaussian type as shown in Fig. |5l The SPM, therefore, can be applied 
and shows us the position of the peak of each term of the sum above; the 
first peak is given by 

-x/v + T-t = (40) 
and the subsequent peaks are given by 

-x/v + 2n{L-a)/v + {2n + l)T-t = (n = l,2, ■■■). (41) 

In particular, the departure times from x = of those peaks are given by 

t = 2n{L-a)/v + {2n + l)T (n = 0, 1, 2, ■ ■ ■)• (42) 

Consecutive departure times are separated by 2(L — a)/i' + 2r. 

A similar discussion for the behavior of the wave fl5B]) in the region a < 
X < L is possible: The coefficients a{E), I3{E) are 

(cosh xa + 2A_ sinh xO')e~^''^^~"'^ 
(cosh xa + iA_ sinh x^f)^ Q-ik(L-a) _|_ sinh^ X'^s*'^^^"'^) 

2n 



i{2nk{L-a)+{2n+l)<l>} 



n=0 
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Figure 6: Motion of the peaks 

— z A+ sinh xae^''^^~"'^ 
(cosh xa + iA_ sinh x^)^ Q-ik{L-a) _|_ sinh^ x'^e*'^^^"'^) 

. , \ 2n+l 
+ SmJl \ gj{2(n+l)fc(L-a)+2{ri+l)(/.} 

r J ' 

which indicate that the multiple peaks of right— and left— going waves 
emerge in the region a < x < L. The right— going peaks depart from 
a; = a at 

t = 2n{L-a)/v + {2n + l)T (n = 0, 1, 2, ■ ■ ■) (43) 
and the left— going peaks arrive at x = a at 

t = 2{n+l){L-a)/v + 2{n + l)T (n = 0, 1, 2, ■ ■ ■)• (44) 

The motions of the peaks thus obtained can be represented as in Fig. 
El from which we can understand that the multiple reflections in the region 
a < X < L cause the existence of the multiple peaks in the transmitted wave 
in the region x > L + a. 

De Leo and Rotelli [7] state that the generalized Hartman effect repre- 
sents an example of an ambiguity in the use of the SPM. They state that 
although the SPM concludes the existence of a generalized Hartman effect 



00 ^ / 

n=0 ' \ 
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if the tranmitted wave is assumed to have a single peak, the same method 
concludes an alternative result if the transmitted wave is assumed to have 
multiple peaks. Admitting the posibility of multiple outgoing peaks, they 
demonstrate the same formula as Eq. fl35l) . 

Following the results obtained in this paper, we have to emphasize two 
points: (a) The generalized Hartman effect is not an example of an ambiguity 
in the use of the SPM, but rather an example of a mathematically incorrect 
application. There is no ambiguity in the use of the SPM. (b) The fact 
that the transmitted wave has multiple peaks is not an assumption. It is an 
inevitable result of a mathematically correct application of the SPM. 

Finally, we have to turn our attention to a report written by Longhi, 
Laporta, Belmonte and Recami f3] in which tunneling optical pulses at 1.5 
fim wavelength through double-barrier periodic fiber Bragg gratings is experi- 
mentally investigated. They claim that the transit time is paradoxically short 
and almost independent of the barrier distance. They interpret their result 
as providing, in the optical context, experimental evidence of the analogous 
phenomenon in quantum mechanics of nonresonant superluminal tunneling 
of particles across two successive potential barriers. 

As we have discussed above, the cautious application of the SPM shows 
that the transmitted wave has an infinite number of peaks and that the time 
each peak emerges at x = L + a depends on the distance L — a. Therefore, 
if the interpretation of the authors of [3] is accepted, we are faced with a 
serious contradiction between theory and physical reality. Resolution of this 
contradiction is very important in order to encourage steady progress in this 
research area. 
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